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[57] ABSTRACT 

A system and method for decoding by identification of the 
most likely phase coded signal corresponding to received 
data. The present invention has particular application to 
communication with signals that experience spurious ran- 
dom phase perturbations. The generalized estimator- correla- 
tor uses a maximum a posteriori (MAP) estimator to gen- 
erate phase estimates for correlation with incoming data 
samples and for correlation with mean phases indicative of 
unique hypothesized signals. The result is a MAP likelihood 
statistic for each hypothesized transmission, wherein the 
highest value statistic identifies the transmitted signal. 


20 Claims, 3 Drawing Sheets 
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5,303,269. 


BACKGROUND OF THE INVENTION 

25 

1. Field of the Invention 

The present invention relates to communications systems 
and, more particularly, to digital signal decoders. 

2. Description of the Related Technology 

30 

Most electronic communication systems in use today 
include a transmitter to transmit an electromagnetic signal 
and a receiver to receive the transmitted signal. The trans- 
mitted signal is typically corrupted by noise and, therefore, 
the receiver must operate with received data that reflects the 
combination of the transmitted signal and noise. Thus, the 35 
receiver receives data y(t) at a time t, where y(t)=s(t)+n(t), 
the sum of the transmitted signal and additive noise. The 
received data equation can be expanded as follows: 

y(/) = A(/) cos[gv + 0(0] + n ( t ) 40 

where A(t) is the signal amplitude, co 0 is the carrier or 
reference frequency, 0 (t) is the time-varying phase function 
and n (t) is noise. 45 

Many of these communication systems require that the 
receiver demodulate information in the received signal 
which depends on proper demodulation of the signal phase 
angle at all times during transmission. The demodulation of 
the signal phase angle is problematic in view of the perva- 50 
siveness of noise. Therefore, for this class of receivers it is 
desirable to optimize phase demodulation, which is equiva- 
lent to optimizing an estimation of the phase function 0 (t). 

Digital communication involves modulation that changes 
discontinuously in accordance with a transmitted coded 55 
signal, e.g., discontinuous phase shifts of a sequence of 
short- duration sinusoids (chip waveforms), as in quadrature 
phase shift keying (QPSK). Each coded signal represents a 
transmitted quantization level or symbol. One example type 
of coded signal is a maximum length binary sequence (W. 60 
W. Peterson, Error-Correcting Codes , MIT Press, Cam- 
bridge, Mass., 1961). Such a coded signal is a binary stream 
which is generated by a decorrelated cyclic shift, i.e., each 
coded signal is designed to be uncorrelated with any another. 
Coded signal transmission is used today in many commu- 65 
nications domains including, for example, digital cellular 
telephone systems. 


Uncompensated multipath, synchronization errors, and 
phase disparity between transmitter and receiver cause a 
received code to have phase values that are different from 
those of the transmitted code, even without additive noise. 
The differences between received and transmitted phase 
samples can be modeled as random perturbations that are 
correlated from sample to sample. A conventional maximum 
likelihood (correlation) receiver does not account for such 
noise-free randomness; received noise-free samples are 
assumed to have the same phase values as the transmission 
for correlation processing. Thus, a need exists for an opti- 
mum receiver/demodulator for digital transmission in the 
presence of such errors. 

The maximum a posteriori (MAP) estimator uses each 
phase value in a hypothesized coded signal as an element of 
the prior mean phase vector 0 m . For the same data, different 
MAP phase estimates are obtained for different hypoth- 
esized signals. The set of MAP phase estimates correspond- 
ing to a given hypothesized coded signal determines a 
reference signal for correlation with the original data. 

The resulting receiver performs better than a conventional 
maximum likelihood processor for phase perturbed signals, 
but a different 0 m vector and associated MAP phase estimate 
is needed for each coded signal hypothesis. For many coded 
signal hypotheses, receiver complexity for MAP estimator- 
correlator processing seems to be prohibitive. Thus, a need 
exists for a simple implementation of a MAP estimator- 
correlator receiver for digital communications. 

SUMMARY OF THE INVENTION 

The problem of optimal demodulation can be addressed 
with a maximum a posteriori (MAP) formulation, with a 
prior phase distribution that has nonzero covariance. Assum- 
ing that phase perturbations have zero mean value, the prior 
phase distribution has a mean phase vector corresponding to 
the transmitted coded signal and a covariance matrix corre- 
sponding to the covariance of the phase perturbations caused 
by multipath, synchronization, and/or phase disparity 
between transmitter and receiver. The prior mean phase 
values in a MAP demodulator thus correspond to the phase 
values in a hypothesized coded signal. 

A correlation receiver that compares received phase shifts 
with hypothesized coded signals (or performs equivalent 
operations on in-phase and quadrature (I,Q) components) is 
suboptimum for random phase perturbations that are not 
caused by additive noise. For noise-independent random 
phase perturbations, the correlator should be replaced by a 
MAP version of a generalized estimator-correlator that first 
forms MAP phase estimates 0 (using received phase shifts 
and hypothesized coded signals) and then correlates these 
estimates with the received phase shifts and hypothesized 
coded signals using the MAP likelihood statistic 

hiApik ) = 

1 N 

= 1— £ {A 2 (iA) + lr(iA)t 2 - 

2a„ 2 i= 1 

2A(i'A)lr(iA)lcos[0jt(iA) - a(i‘A)]} - 
i N N 

4 - s X [ 0 fcO‘A) - 

z i=l j= 1 

where p (rl0 ft ) and p (0 fc l0 mfc ) are conditional probability 
density functions, r(iA)=lr(iA)!exp[ja(iA)], i=l, . . . , N is the 
sampled received signal, C 0 is the phase covariance matrix, 
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the index k refers to the k th hypothesized coded signal, 
0 mJfc is the mean phase vector for the k //! coded signal and the 
MAP estimates 

0 fc =[6 A (A) . . . 0*(NA)] 7 

satisfy the equation 

0*(/A) - 6ffljt(iA) = g~ 2 A 2 lr(/A)l/?ei(r -j)A\sin\a(JA) - 0jt(/A)], 

7=1 

i=l, . . . , N, where R G [(i-j)A] corresponds to the element in 
row i and column j of the phase covariance matrix C 0 . When 
R 0 [(i-j)A] equals zero, the MAP estimator yields the prior 
mean phase estimates 0 m and the generalized estimator- 
correlator becomes an ordinary correlator. 

A simpler estimator-correlator configuration was intro- 
duced in the original patent application in conjunction with 
a classifier that discriminates between signals with different 
phase covariance functions, R 0 (.). The present invention 
discriminates between phase-perturbed coded signals with 
different mean phase values corresponding to the transmitted 
phase values of each signal. The coded signal discriminator 
uses a generalized version of the estimator-correlator con- 
figuration, as described in the above equations. 

The complexity of a generalized estimator-correlator 
receiver can be reduced by exploiting the cyclic property of 
many block codes, e.g., Reed-Solomon codes. If a length N 
code is cyclic, a sequence of N-l cyclic shifts of a given 
coded signal can generate all the other signals in the code. 
For a unit cyclic shift to the right, the last element in a shift 
register is moved to the first position, and all other elements 
are shifted to the right by one position. 

A single estimator-correlator can test all N coded signals 
in a length N cyclic code where the convergence time of the 
MAP estimator is assumed to be less than or equal to the 
sampling interval A. This requirement is also implicit in a 
block processor that moves a block by one sampling interval 
for each MAP phase estimate. After NA seconds, the N 
sample input buffer is filled with samples corresponding to 
the next coded signal. During the same time, N cyclic shifts 
of a hypothesized coded signal yield N sets of MAP phase 
estimates to be used in a generalized estimator-correlator. 
The largest estimator-correlator output designates the best 
hypothesis for the current coded signal. Phase samples of the 
next received signal are transferred to the estimator-correla- 
tor via a parallel set of connections between each element of 
the input buffer and the corresponding element of the data 
shift register. A conventional maximum likelihood demodu- 
lator can exploit the same principle as the MAP demodula- 
tor. 

In summary, the present invention utilizes phase values of 
a hypothesized coded signal as the prior mean phase vector 
in a generalized MAP estimator-correlator demodulator for 
digital communication with random phase perturbations. 
Random phase perturbations are introduced by synchroni- 
zation (timing) errors, uncompensated multipath, and 
uncompensated phase shift between transmitter and receiver. 
The present invention also utilizes recirculating shift regis- 
ters for efficient implementation of ML and MAP demodu- 
lators for cyclic codes. 

These and other objects and features of the present 
invention will become more fully apparent from the follow- 
ing description and appended claims taken in conjunction 
with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of one presently preferred 
embodiment of a maximum a posteriori (MAP) phase angle 
5 demodulator using a Hopfield network. 

FIG. 2 is a block diagram of a MAP decoder of the present 
invention, including a MAP phase estimator as shown in 
FIG. 1, for a length N cyclic code, assuming uniform 
amplitude |rl=[I 2 +Q 2 ] 1/2 

10 FIG. 3 is a block diagram of a set of generalized cross 
correlators for MAP decoding, i.e., for determining the most 
likely transmission for K coded-signal hypotheses. 

15 DETAILED DESCRIPTION OF THE 

PREFERRED EMBODIMENTS 

FIG. 1 illustrates one preferred embodiment of a maxi- 
mum a posteriori (MAP) demodulator 100. The background 
of a general implementation for the MAP demodulator 100 
20 can be found in the co-pending parent application, now U.S. 
Pat. No. 5,303,269, which is hereby incorporated by refer- 
ence. 

Essentially, the MAP demodulator 100 receives a phase 
25 modulated signal such as, for example, an FM broadcast 
signal. A sampling circuit 102 periodically samples the 
integrated products of the received signal y(t) with two local 
oscillator signals, cos(co 0 t) and sin(eo 0 t), where co 0 is the 
carrier frequency. The products of y(t) cos oo 0 and y (t) sin 
3Q co 0 are integrated over an interval A, and the integrator 
outputs are sampled at intervals of jA, j=l, . . . , N, to 
generate sampled in-phase and quadrature components, I(jA) 
and Q(jA), j=l, . . . , N. The sampled components are fed into 
an arctangent circuit 104 and an amplitude circuit 106. The 
35 arctangent circuit 104 performs the function arctan(Q/I), 
where the received signal comprises in-phase (I) and quadra- 
ture (Q) components, to provide a time series of phase 
samples a(jA). The amplitude circuit 106 outputs a time 
series of sample amplitudes ly(jA)i. The N sample phases 
4Q and amplitudes from the circuits 104,106 are fed into a 
4 Hopfield MAP phase estimator 108. The MAP estimator 108 
also receives N predicted mean phase values 0 m (jA), j=l, . 

. . , N. The MAP estimator 108 provides a block of N 
estimated phase values 0(kA), k=l, . . . , N at its output. One 
45 of the N samples in the output block is shown at the output 
of the phase estimator 108. Thus, the desired goal of the 
MAP demodulator 100 is to achieve an accurate estimate of 
the phase sequence of the received signal y(t). 

A coded signal hypothesis generator 110 generates the 
50 prior mean phase for a particular coded signal k, k=l . . . K, 
which is fed into the MAP phase estimator 108. The gen- 
erator 110 may also generate code word amplitudes. How- 
ever, for the purpose of explaining the general demodulator 
100, this aspect is not shown in FIG. 1 as it will be discussed 
55 below with respect to FIG. 3. 

The presently preferred MAP estimator 108 is structurally 
arranged as follows. A summing amplifier 112 accepts the 
phase samples from the arctangent circuit 104 and subtracts 
each sample from a corresponding predicted mean phase and 
60 provides the result to a Hopfield network 114. Phase differ- 
ences {a(jA)-0 mA .(jA)} are constrained to lie between -n and 
n by addition of 2 n, 0, or -2k to relevant a-samples. The 
Hopfield network 114 also receives amplitude samples cor- 
responding to the input signal after, optionally, processing 
65 the signal via a scaling circuit 116. Shift register 118 and the 
multipliers 120 which multiply signal amplitude samples by 
corresponding samples of the autocovariance function 
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R 0 [(k-j)A]. The output of the Hopfield network 114 is fed 
into a summing amplifier 122 to result in the phase estimate 
0(kA), the end product of the demodulator 100 . 

For analog communication, maximum a posteriori (MAP) 
phase estimation is sufficient for demodulation. Digital 5 
systems, however, require a hypothesis test to decide 
between transmitted coded signals, bits, or quantized chip 
waveform phase shifts. MAP phase estimates can be used to 
replace a computationally expensive optimum hypothesis 
test that involves multidimensional integration of the prod- 10 
uct of two multivariate probability distributions. The gen- 
eralized estimator-correlator requires that the vector 0* of 
MAP phase estimates be correlated with both the prelimi- 
nary (maximum likelihood) phase estimate vector a (as in an 
estimator-correlator) and with the prior or reference phase 
vector Q mk 

To understand the structure and function of the present 
invention it is helpful to see the mathematical derivation of 
a coded signal hypothesis test. The hypothesis test is imple- 
mented in order to decide which of K transmitted phase 
vectors {0 mJfc , k=l, . . . , K} is most likely, based on a 20 
received complex data vector E and a prior phase distribu- 
tion describing phase uncertainty that could be caused by a 
rapidly time varying, unequalized communication channel 
and/or system errors (e.g., imperfect synchronization and 
spurious phase shifts). If all the coded signal mean phase 
vectors {0 m J are equally likely a priori, the most likely 
0 mJt maximizes the likelihood function p(rl0 mJt ), k=l, . . . , K. 

The hypothesis test accounts for channel/system induced 
uncertainty in the received 0-vector by integrating over all 
possible received 0-vectors, weighted by their prior prob- 30 
abilities; 

P0%,*)=J. . . Ip(rifi) pffliejde. (1) 

The multidimensional integral in Equation (1) integrates 35 
each component of 0 over a 271 interval, where 

0 r =[0(A)0(2A) . . . 0(NA)] (2) 

and the sampling interval A represents the duration of a 
constant frequency, phase shifted “chip” waveform in a PSK 40 
signal. 

The multivariate probability distributions in Equations 

(1) are 

P (rl 0)=( (2 tc) /W 2 I C„ 1 1/2 ]~ 1 exp l-{V 2 )(r-s 0 ) T C n -\r- SQ )\ (3) 45 
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A simpler implementation of a hypothesis test based on 
the likelihood function in Equation (1) is based on the 
following observation: In additive, white noise, the 
Q mk - vector that maximizes the likelihood function p(rl0 wfc ) in 
Equation (1) is also the 0^- vector that maximizes p(rl 
where 0* is the MAP 0-estimate corresponding 

to Qmk- ~ 

Let phase sequence Q mj be the most likely hypothesis 
obtained via multivariate integration; 


J...Jp(H0)p(eie m ;)<ie = 


max 




where C n ~o n 2 l and the integration limits are 


(7) 


{ a(rA)-7i, A(i A)+je; i= 1 , . . . , N}. (8) 

Let 0 m „ be the chosen hypothesis obtained by substituting 
the MAP phase estimate into the integrand in Equation (7); 


p(M n ) P (m m )\ = 


max 

®mk 


] (pye t )p(e*i0 mt )i 


(9) 


v / * 

where the MAP phase estimate 0* is the 0 mr vector that 
maximizes p(rl0)p(0l0 TO J. It can be shown that 0 mn equals 
Q mJ , i.e., that the hypothesis based on the MAP estimate is 
the same as the most likely hypothesis. 

Replacement of p(rI0 wJk )=J. . . Jp(rI0) p(0l0 m *)d0 in Equa- 
tion (7) by the simpler statistic in Equation (9) is important 
because evaluation of the multivariate distributions for all 
0- samples and integration over all such samples is compu- 
tationally burdensome unless C 9 is a diagonal matrix. Deter- 
mination of the MAP phase estimate 0* for each transmitted 
phase hypothesis 0 mA: may also be infeasible unless an 
analog Hopfield network (or a parallel set of K such net- 
works for 0 m/ , . . . , 0 mAr ) is used to solve the equations 


, N ~ (10) 

e*(rA) - 0„*(/A) = a n 2 Z A(/A)lr(/A)l/?ol(*' - ;)A]sin[a(/A) - 0*(/A)], 

H 

i - 1, . . . , N, k = 1, . . . , K. 

Computation of the MAP based test statistic p(rl0*)p(0*l 
0 mJt ) involves a generalized correlation operation. Maximi- 
zation of the statistic over {0 m *; k=l, . . . , K} is equivalent 
to finding the 0 mJt vector that maximizes 


lnfp(rie Jfc )p(0 jfc l0 mJfc )l = 


(ID 


and 


p(Sia„ A )=[(27c)^lC o l 1/2 ]- 1 cxp[-(^)(e-0 mJfc ) :r Ce- 1 CO-0 W i fc )] (4) 

where the \ th components of, the sampled received data time 
series, and s 0 , the transmitted signal function, are 


l — Z {A 2 (iA) + lr(/A)l 2 - 

2ct„ 2 i= 1 

2AIKrA)lcos[0 k (/A) - a(/A)]} - 

-J- .1 .1 [0*(i4)-e^(;A)]|Q'ye*(;A)-eM(/'A)]. 

z i-l 7=1 


r(rA)=lr(rA)lexp[za(*A)] (5) 

and 

So(i‘A)=A([A)exp[j0(*'A)] (6) 

respectively. 

A maximum likelihood phase estimate computes the 
phase of s0 in Equation (6) that maximizes the conditional 
probability p(ri0) of signal r given phase 0 in Equation (3), 
and this phase vector is the same as a in Equation (5). The 
phase covariance matrix C 0 in Equation (4) describes the 
deviation of received phase values from their expected 
values, and is associated with unequalized or randomly time 
varying channel effects or system errors. 


The first expression on the right hand side of Equation (11) 
55 is maximized by maximizing the output of an estimator- 
correlator that computes 

Re | S \r(iA)ti a MxA(iA)e-flkm J (U) 

60 The second expression on the right hand side of Equation 
(11) is proportional to a different correlation operation 
involving the MAP phase estimate, i.e., 

8jfc 7 Ce ’ll**- (13) 

65 

The estimator-correlator in the original patent application 
can be distinguished from that of the present invention, since 
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it discriminates between processes with different covariance 
functions, e.g., speech vs. music, rather than between 
sequences with different prior mean phase vectors, e.g., 
different coded signals, as considered here. 

FIG. 2 illustrates a block diagram of a MAP decoder 128 , 5 
of the present invention, for a length N cyclic code, assum- 
ing uniform amplitude lrl=[I 2 +Q 2 ] 1/2 . The coded signal deci- 
sion corresponds to the shift in the recirculating signal shift 
registers that yields maximum correlator response \ MAP , 
which is the MAP likelihood function. 10 

The MAP decoder 128 shown in FIG. 2 is an efficient 
design for use in conjunction with cyclic codes. The MAP 
decoder for cyclic codes 128 is structurally connected as 
follows. A sampled received or data signal having a 
sequence of digitally coded phase shifts is fed into an 15 
arctangent circuit 130 , which provides sampled phase input 
to a data sample input buffer 132 . Note that for the illustrated 
MAP decoder 128 , the hypothesized amplitude A and 
sampled received signal amplitude lr(iA, I are assumed to be 
constant and are, therefore, not shown. 20 

The input data from the buffer 132 is fed into a data 
sample shift register 134 and the MAP phase estimator 108 . 
The data shift register 134 feeds sample data to a generalized 
cross correlator 136 . The generalized cross correlator 136 
also receives MAP phase estimates 0* from a MAP phase 25 
shift register 138 , receivably connected to the output of the 
MAP phase estimator 108 , and prior mean phases Q MK , 
representative of the k possible coded signals, from a 
recirculating prior mean phase shift register 140 . The output 
of the generalized cross correlator 136 , is the MAP likeli- 30 
hood function result correlating the input phase sequence 
with each of the K hypothesized coded signals. Although not 
shown in FIG. 2, ultimately the MAP decoder 128 must 
compare the K correlations and identify the input phase 
sequence to be the coded signal with the highest value 35 
statistic. 

The MAP estimator 108 can be applied to demodulation 
of digital communication signals if the block of phase 
samples coincides with a coded signal, i.e., a sequence of 
phase shifts corresponding to a transmitted symbol or digit. 40 
The prior mean phase values 0 m correspond to the phase 
samples of the hypothesized coded signal. The autocovari- 
ance function R 0 [(i-j)A] corresponds to the covariance 
matrix C e of the prior phase distribution. 

The MAP demodulator is a generalized version of the 45 
estimator-correlator (H. L. Van Trees, Detection , Estimation, 
and Modulation Theory, Part /, Wiley, New York, 1968) as 
in Equation (11) above, such that MAP phase estimates 
{0(kA)| are correlated not only with corresponding prelimi- 
nary phase estimated {ce(kA)} extracted from incoming data 50 
as in Equation (12) but also with C o _1 0,„. The hypothesized 
coded signal with the largest generalized estimator-correla- 
tor output 1 MAP is chosen as the most likely transmission. 

The MAP demodulator 100 (FIG. 1 ) is a generalization of 
a standard maximum likelihood (ML) demodulator. The ML 55 
demodulator correlates estimated phase values (or the cor- 
responding in-phase and quadrature components I,Q) with 
hypothesized phase values 0 m (or the corresponding I,Q 
components) for each hypothesized coded signal. The 
hypothesized coded signal with maximum correlator output 60 
is chosen as the most likely transmission. MAP and ML 
demodulators are identical in the limit as R 0 [(H)A] 
becomes very small for all i-j, corresponding to a prior 
phase distribution that is concentrated at 0=0 m . In this case, 
the reference signal becomes known instead of random, 65 
leading to a standard correlationjeceiver. The corresponding 
MAP phase estimate becomes 0=0 m . 
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An ML receiver, which is not shown in FIG. 2, would 
require a cross correlator between shift register 140 and shift 
register 134 . This correlation actually occurs in FIG. 2 when 
phase uncertainty as represented by C Q becomes very small. 

MAP demodulation is superior to ML when the transmit- 
ted phase is randomized, even when no noise is added. 
Randomized transmitted phase shifts may have correlated, 
random phase perturbations. Such perturbations are intro- 
duced by synchronization (timing) error, an uncompensated 
multipath component, or random phase errors introduced by 
local oscillator phase drift or transmitter/receiver motion 
(Doppler) effects. 

The MAP estimator 108 for M-ary PSK (phase-shift 
keyed signals with M different permissible chip waveform 
phase shifts) uses each phase value in a hypothesized coded 
signal as an element of the prior mean phase vector 0 m . The 
estimator also uses corresponding amplitude and phase 
measurements obtained from received data. For the same 
data, different MAP phase estimates are obtained for differ- 
ent hypothesized coded signals. Each of these estimates is 
substituted into a reference signal for correlation with the 
original data. The resulting receiver performs better than a 
conventional ML processor for phase perturbed signals, but 
a different 0 m vector and associated MAP phase estimate is 
needed for each coded signal hypothesis. If many coded 
signal hypotheses are implemented in parallel, receiver 
complexity for MAP estimator-correlator processing seems 
to be prohibitive. 

A solution to the complexity problem is obtained by 
exploiting the cyclic property of many relevant block codes, 
including M-ary Reed-Solomon codes (J. G. Proakis, Digital 
Communications , 2nd ed., McGraw-Hill, New York, 1989, 
pp. 429-432; W. W. Peterson, Error-Correcting Codes , MIT 
Press, Cambridge, Mass., 1961, pp. 137-182). If a length N 
code is cyclic, a sequence of N-l cyclic shifts of a given 
coded signal (e.g., s } ) can be used to generate all the other 
signals (s 2 , . . . , s N ) in the code. For a unit cyclic shift to the 
right, the last element in a shift register containing S 2 is 
moved to the first position, and all other elements are shifted 
to the right by one position. The length of each coded signal 
thus corresponds to the number of different codes N. 

Any coded signal in a cyclic code of length N can be 
decoded with a single ML correlator by using the following 
procedure: 

(1) Obtain the in-phase and quadrature components of 
each sample of the received signal, construct the correspond- 
ing amplitude and phase variables lr(jA)KI 2 +Q 2 ) 1/2 and 
oc(jA)=tan -1 (Q/I), and store these quantities in an input 
buffer. After samples corresponding to an unknown coded 
signal of length N have been obtained, transfer the samples 
to a data shift register. The data shift register is to be used 
for correlation of stored data samples with corresponding 
signal samples that are stored in a reference signal shift 
register. 

(2) Assuming that the first coded signal corresponds to a 
sampled signal of the form s 1 (iA)=Aexp[j0 m (iA)], i=l,2. . . 
, N, store the 0 m vector with elements 0 m (A), . . . , 0 m (NA) 
in an N-sample recirculating signal shift register. 

(3) Cross correlate samples in the data shift register with 
signal samples obtained by cyclic shifting of 0 m , for shifts of 
0, . . . , N-l. For k shifts, the cross correlation function is 
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£ lr(iA)ld“( ,A Me~ j0mA(,Ai 

h«,(fc) = 02 

N N 1 

£ I r(z A)l 2 £ [A I 2 
i=l i=l J 

where 8 m *. is a version of 8 m that has been cyclically shifted 
k times. 

(4) Determine the shift m that yields maximum correla- 
tion, 

Re{lMi{m)} = max 
0*k*N-l 

and output the corresponding decision, viz., that coded 
signal m+1 was transmitted. 

As shown in FIG. 2, the MAP estimator- correlator 128 
uses samples lr(iA)l and a(iA), i=l, . . . , N, in the data shift 
register 134 and samples 8 wJfe (iA), i=l, . . . , N in the signal 
shift register 140 to compute via the MAP estimator 108 a 
vector of MAP phase estimates 9 fc corresponding to hypoth- 
esized coded signal k+1. The MAP phase estimates are 
obtained via steepest descent solution of a nonlinear integral 
equation with a Hopfield network that has convergence time 
t c . The MAP phase estimates 8* are stored in a reference 
signal shift register 138 for a generalized correlator that 
computes the MAP likelihood function for the k„, coded 
signal as follows: 

W*Hnlp(rIg*)p (§*10^)1 04) 

where the right hand side of Equation (14) is the same as 
Equation (11). 

The MAP decoder 128 shown in FIG. 2 will process data 
in real time if N MAP phase vector estimates can be obtained 35 
within one coded signal duration, i.e., if Nt c NA, or t c ^A, 
where t c is the convergence time of the Hopfield MAP 
estimator. A delay of Nt c seconds is needed to allow the 
recirculating shift register 140 to cycle through all possible 
cyclic shifts. If t c ^A, this delay is less than or equal to the 40 
time NA required to fill the N-sample input buffer with data 
samples pertaining to the next transmitted signal. 

If a VLSI Hopfield network circa 1989 (A. Moopenn and 
A. P. Thakoor, “Programmable synaptic devices for elec- 
tronic neural nets,” Proc. 5th IASTED Int’l Conf. on Expert 45 
Systems and Neural Networks, Honolulu, Hi., Aug. 1989, 
pp. 36-40) with convergence time between one and ten 
microseconds is used, then t c ^A if maximum channel 
bandwidth B=1/(2A) is between 50 kHz and 500 kHz. 
Processing speed and corresponding channel bandwidth can 50 
be increased by using newer, faster Hopfield implementa- 
tions or by processing a sequence of N-sample data blocks 
in parallel. For sufficiently short convergence time, the MAP 
decoder 128 shown in FIG. 3 requires only one Hopfield 
network and correlator for MAP demodulation of a code of 55 
arbitrary size. The only disadvantage is a processing delay of 
Nt c seconds, which increases linearly with the code length 
N. 

FIG. 3 is a detailed block diagram of the generalized 
estimator-correlator 128, of the present invention, particu- 60 
larly showing the structure of the generalized cross correla- 
tor 136. The generalized cross correlator implements the 
statistic in Equation (11). Note that the construction of the 
decoder 128 requires K parallel circuits. Thus, there is one 
circuit per hypothesized coded signal. Note that for the 65 
special case of cyclic codes, the decoder 128 would only 
require one such circuit, as shown in FIG. 2. 
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In FIG. 3, in-phase and quadrature components I(iA) and 
Q(iA) are extracted from the received signal, and the result- 
ing I,Qlr(iA)i=[I 2 +Q 2 ] 1/2 and a(iA)=arctan(Q/I) values may 
be stored at the data sample input buffer 132. The sampled 
5 signal is then distributed to the K parallel circuits at the data 
sample correlators 152 and the MAP phase estimators 108. 
(The first, second and kth circuits are indicated in FIG. 3 by 
the use of a, b, and c suffixes, respectively.) 

Each coded signal hypothesis comprises an amplitude A 
10 and a mean phase 9 W which are stored in a hypothesized 
coded signal data register 154. The coded signal data register 
154 feeds the data sample correlator 152, the MAP phase 
estimator 108 and a mean phase correlator 156. The MAP 
15 phase estimator 108 provides a block of phase estimates 9 to 
the correlators 152, 156. 

The correlated values from the correlators 152, 156 are 
fed into a summing amplifier 158. The result of each 
summation is the MAP likelihood statistic 1 MAF . Each of the 
20 K calculated statistics is fed into a comparator where the 
highest value identifies the received data as one of the 
hypothesized coded signals. 

In the generalized estimator- correlator 128 of FIG. 3, the 
sampled data at time iA is y(iA)=lr(iA)lexp[ja(iA)]. I(iA) and 
25 Q(iA) are extracted as previously described. 

Ir(/A)l 2 = P(iA) + Q 2 (iA) a(/A) = tairhgO'Ay/O'A)] 

cos[a(zA)] = I(iA)/\r(iA)\ sin[a(/A)] = g(*A)/lr(iA)]. 

30 The Hopfield MAP phase estimator 108 of FIG. 1 has 
Re[(H)A]=[C0](/- ReKH) A] is the element in row i and 
column j of the phase covariance matrix C Q , where C e =E{( 
0-0 m )(0-8 w ) 7 }. The simplest technique for estimating R e (.) 
or the phase covariance matrix is via decision feedback. If 
the column vector 8^ is the most likely coded signal and if 
9* is the corresponding MAP phase vector, then an estimate 
of C 0 is the sample covariance matrix^ corresponding to a 
weighted average of past versions of (8-9 mA: )(0-9 mA: ) :r . The 
output of the Hopfield MAP phase estimator is a vector 
8 fc of MAP phase estimates corresponding to the k /ft coded 
signal hypothesis, i. e. , to the symbol or coded signal A 
*(iA)exp [j8 mJfc (iA) ], i=l, . . . , N. The a, 8 correlator 152 
computes 

N ( 15 ) 

- £ HrOA)!exp[/'a(iA)] - 
i'=l 

N 

A(iA)exp[/0*(/A)]l 2 = £ {A(zA)[/(zA)cos0*(*A) + 
r=l 

e(*A)sine k (/A)] - lr(zA)l 2 - A 2 (rA)} 

and the 8,8 m correlator 156 computes 
-(i-an^Ce- 1 ^ 

( 16 ) 

It will be recognized that the present invention can be 
implemented in a variety of different ways. As but one 
example, the demodulator could be implemented using 
integrated circuit components, including a significant 
amount of digital logic, but the Hopfield network would 
preferably be constructed using analog circuitry, requiring 
appropriate D/A and A/D conversion, to achieve a fast rate 
of convergence of course, standard engineering techniques 
allow for implementations tending to emphasize computer, 
optical, electrical, or other methods, with engineering and 
economic trade-offs restricting each implementation as a 
practical matter. 



While the above detailed description has shown, 
described and pointed out the fundamental novel features of 
the invention as applied to various embodiments, it will be 
understood that various omissions and substitutions and 
changes in the form and details of the device illustrated may 5 
be made by those skilled in the art, without departing from 
the spirit of the invention. 

What is claimed is: 

1. A maximum a posteriori decoder, comprising: 

a data sample input buffer receiving a plurality of samples 10 
of a received signal; 

a data sample shift register receivably connected to the 
data sample input buffer; 

a mean phase shift register storing a plurality of mean 
phases each indicative of a unique coded signal; 15 

a maximum a posteriori phase estimator receivably con- 
nected to the data sample input buffer and the mean 
phase shift register, said estimator generating a plural- 
ity of phase estimates; 

a phase estimate shift register receivably connected to the 
phase estimator; and 

a generalized cross correlator receivably connected to the 
data sample shift register, the mean phase shift register 
and the phase estimate shift register, said cross corr- 25 
elator generating a likelihood statistic for each unique 
coded signal. 

2. The maximum a posteriori decoder defined in claim 1, 
wherein each of the samples comprises the phase of a chip, 
where a chip is defined as a short, constant frequency 30 
waveform comprising part of a received signal. 

3. The maximum a posteriori decoder defined in claim 1, 
wherein each of the samples comprises an amplitude lr(iA)l 
and a phase a(iA), i=l, . . . , N. 

4. The maximum a posteriori decoder defined in claim 1, 35 
wherein the mean phase shift register cycles through all 
possible cyclic shifts in Nt c seconds, wherein N is the 
number of mean phases and t c is the convergence time of the 
maximum a posteriori phase estimator. 

5. The maximum a posteriori decoder defined in claim 1, 40 
wherein the generalized cross correlator computes the maxi- 
mum a posteriori likelihood function for the k //2 coded signal 

as follows: 

where 

^the complex data vector associated with the received 
signal; 

G^^the mean phase vector for the k th coded signal; 50 

G^the vector of chip phase shifts corresponding to a k" 1 
hypothesized transmitted signal for direct sequence, 
spread spectrum communication. 

6. The maximum a posteriori decoder defined in claim 1, 

wherein the generalized cross correlator includes: a data 55 
sample correlator for correlating each phase sample and 
phase estimate, a mean phase correlator for correlating each 
mean phase and phase estimate, and a summer receivably 
connected to the data sample correlator and the mean phase 
correlator to output the likelihood statistic. 60 

7. The maximum a posteriori decoder defined in claim 1, 
wherein the maximum a posteriori phase estimator includes 
a Hopfield network. 

8. A method of maximum a posteriori (MAP) decoding of 

a signal to generate a maximum a posteriori likelihood 65 
function for a plurality of coded signals, the method com- 
prising the steps of: 


receiving a signal; 

sampling in-phase and quadrature components of the 
received signal; 

estimating MAP phase estimates 0 of the received signal 
from the sampled signal components; 
obtaining prior mean phases B mk for a k! h coded hypoth- 
esized transmitted signal; and 
cross-correlating 0 and B mk to obtain a MAP likelihood 
statistics for each k /;i coded signal. 

9. The method of MAP decoding defined in claim 8, 
wherein each MAP phase estimator is parameterized by a 
time series of phase samples a(jA), and a time series of 
sample amplitudes ly(jA)l, j=l, . . . , N. 

10. The method of MAP decoding defined in claim 8, 
wherein the cross-correlating step includes the steps of: 

correlating each data sample with a corresponding one of 
the phase estimates; 

correlating each prior mean phase with a corresponding 
one of the phase estimates; and 
summing the data sample and prior mean correlations to 
calculate each MAP likelihood statistic. 

11. The method of MAP decoding defined in claim 10, 
additionally comprising the step of comparing the MAP 
likelihood statistics to identify the highest MAP likelihood 
statistic. 

12. A maximum a posteriori decoding system, compris- 
ing: 

means for receiving a plurality of samples of a received 
signal; 

means for storing the received signal samples receivably 
connected to the receiving means; 
first shift means for storing and shifting a plurality of 
sampled phases receivably connected to the signal 
storage means; 

second shift means for storing and shifting a plurality of 
mean phases each indicative of a predetermined coded 
signal; 

means for generating a plurality of phase estimates receiv- 
ably connected to the signal storage means and the 
mean phase storage means; 

third shift means for storing and shifting a plurality of 
phase estimates receivably connected to the phase 
estimate generating means; and 
means, receivably connected to the phase estimate storage 
means, the sampled phase storage means, and the mean 
phase storage means, for generating a likelihood sta- 
tistic for each predetermined coded signal. 

13. The maximum a posteriori decoding system defined in 
claim 12, additionally comprising means for comparing the 
likelihood statistic for each unique coded signal hypothesis 
and determining the highest likelihood statistic. 

14. A maximum likelihood correlator, comprising: 
a sampling circuit receiving a coded signal; 

an arctangent circuit receivably connected to the sampling 
circuit; 

an amplitude circuit receivably connected to the sampling 
circuit; 

an input buffer receivably connected to the arctangent 
circuit and the amplitude circuit; 
a data shift register receivably connected to the input 
buffer; 

a recirculating signal shift register storing a prior mean 
phase 0 W vector with elements 0 m (A), . . . , 0 m (NA); 
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a cross-correlator receivably connected to the data shift 
register and the recirculating signal shift register, said 
cross-correlator generating a likelihood statistic of the 
coded signal with respect to each prior mean phase. 

15. The maximum likelihood correlator defined in claim 5 
14, wherein the arctangent circuit calculates the phase 
variable a(jA) of the coded signal, j=l, . . . , N. 

16. The maximum likelihood correlator defined in claim 
14, wherein the amplitude circuit calculates the amplitude 
variable fr(jA)l of the coded signal, j=l, . . . , N. 

17. The maximum likelihood correlator defined in claim 
14, wherein the input buffer stores the phase variable a(jA) 
and the amplitude variable lr(jA)l, j=l, . . . , N. 

18. The maximum likelihood correlator defined in claim 
14, wherein the data shift register receives the phase variable 15 
a(jA) and the amplitude variable lr(jA)l from the input buffer 
after N samples of the coded signal have been obtained, j=l, 

. . . , N. 

19. The maximum likelihood correlator defined in claim 
14, wherein the cross correlator correlates samples in the 
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data shift register with signal samples in the recirculating 
signal shift register using the function fork shifts as follows: 

N 

Z \r(iA)\Qi am Ac-j°mk^ 

hulk) — ' YJ2 

N N 1 

Z lr(*A)1 2 Z \A\ 2 

i= 1 i= 1 J 

1Q where 0 mjt is a version of 0 m that has been cyclically shifted 
k times, 0=k^N— 1. 

20. The maximum likelihood correlator defined in claim 
19, wherein the shift m, m=k, that yields maximum corre- 
lation is calculated as follows: 

, Re{l ML {k)}. 

O^kfkN - 1 

* * * * * 



